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Abstract 

Linear  system  identification  of  complex  nonlinear  systems,  such  as  large  space  structures, 
can  be  difficult  because  such  systems  are  often  lightly  damped,  have  dense  modal  spacing, 
respond  to  disturbances  over  a  large  frequency  bandwidth,  and  exhibit  nonlinear  responses. 
Many  excitation  signals  used  to  collect  frequency  response  data  for  linear  system  identifica¬ 
tion  are  poorly  suited  to  systems  that  exhibit  nonlinear  responses.  Specifically,  random  noise, 
burst  random  noise,  pulse-impact,  and  multisine  excitation  signals  can  create  undesired  non¬ 
linear  distortions  in  the  frequency  response  data  that  are  generally  indistinguishable  from 
the  linear  frequency  response.  Thus,  sine  dwell  (also  called  stepped  sine)  excitation  signals 
are  often  used  to  obtain  frequency  response  data  since  sine  dwell  excitation  signals  allow 
the  base  tone  response  to  be  separated  from  higher  order  harmonics.  However,  sine  dwell 
testing  can  be  very  time  consuming.  In  this  paper,  we  present  a  novel  sequence  of  multisine 
excitation  signals  that  are  more  time  efficient  than  sine  dwell  excitation  while  eliminating  the 
effects  of  nonlinearity  induced  harmonics  from  the  frequency  response  data.  The  sequential 
multisine  excitation  signal  is  demonstrated  on  the  Deployable  Optical  Telescope  testbed  at 
Air  Force  Research  Laboratory. 


1  Introduction 


Large  space  structures  often  require  vibration  suppression  to  meet  performance  require¬ 
ments.  Strict  performance  requirements  over  large  frequency  bandwidths  often  prevent  pas¬ 
sive  vibration  isolation.  Active  control  can  be  used  to  achieve  high  levels  of  disturbance 
rejection.  However,  many  active  control  methods  require  high  fidelity  linear  dynamic  models 
that  approximate  the  underlying  dynamics.  In  addition,  high  fidelity  models  are  required  for 
performance  prediction.  Identifying  linear  dynamic  models  of  large  space  structures  can  be 
difficult  because  large  space  structures  are  lightly  damped,  have  dense  modal  spacing,  have 
large  dynamic  range,  respond  to  disturbances  over  a  large  frequency  bandwidth,  and  may 
exhibit  nonlinear  responses. 

There  are  numerous  excitation  signals  used  for  linear  system  identification  in  the  fre¬ 
quency  domain,  including  random  noise,  burst  random  noise,  pulse-impact,  multisine,  and 
sine  dwell  (also  called  stepped  sine)  [1-8].  Many  excitation  signals  used  to  collect  frequency 
response  data  are  poorly  suited  to  systems  that  exhibit  nonlinear  responses.  Specifically,  the 
nonlinear  responses  to  these  excitations  creates  undesired  artifacts  in  the  frequency  response 
data.  Furthermore,  these  undesired  artifacts  are  indistinguishable  from  the  linear  frequency 
response  and  can  significantly  impact  the  characteristics  of  a  linear  model  identified  using 
frequency  domain  system  identification  methods. 

In  particular,  it  is  difficult  to  use  broadband  excitation  signals  to  obtain  linear  frequency 
response  data  for  a  system  that  exhibits  nonlinear  distortions  due  to  harmonics.  In  [4,  5] , 
it  is  suggested  that  a  multisine  excitation  signal  may  be  designed  to  reduced  or  eliminated 
nonlinear  distortions.  Specifically,  nonlinear  distortions  may  be  eliminated  by  removing 
the  appropriate  harmonics  from  a  multisine  excitation  signal.  For  example,  if  there  is  a 
nonlinearity  that  excites  the  odd  harmonics  of  frequency  fi,  then  the  odd  harmonics  of 
/i  should  be  removed  from  the  multisine  excitation  signal.  However,  this  method  requires 
certain  knowledge  of  the  system’s  nonlinearities,  that  is,  we  must  know  the  harmonics  that 
are  being  excited.  In  addition,  the  method  suggested  in  [4, 5]  is  difficult  to  implement  when 
the  system  has  dense  modal  spacing  and/or  there  are  many  base  tones  for  which  higher- 
order  harmonics  are  being  excited.  Therefore,  sine  dwell  excitation  signals  are  often  used  to 
eliminate  nonlinear  distortions.  However,  sine  dwell  testing  can  be  much  less  time  efficient 
than  broadband  excitation  signals  [7].  In  practice,  sine  dwell  testing  is  often  abbreviated, 
resulting  in  lower  quality  data. 

In  the  present  paper,  we  consider  a  novel  sequence  of  multisine  excitation  signals  designed 
to  obtain  system  identification  data  over  a  large  frequency  bandwidth,  while  eliminating 
nonlinearity  induced  harmonics  from  the  frequency  response  data.  Furthermore,  the  method 
presented  in  this  paper  is  more  time  efficient  than  sine  dwell  testing. 

In  Section  2,  we  formulate  the  excitation  signal  design  problem.  The  properties  of  random 
noise,  burst  random  noise,  pulse-impact,  multisine,  and  sine  dwell  excitation  signals  are 
reviewed  in  Sections  3-5.  Section  6  compares  the  test  times  associated  with  multisine  and 
sine  dwell  excitation  signals.  Sequential  multisine  excitation  signals  are  presented  in  Section 
7.  Section  8  provides  a  sequential  multisine  implementation  for  multi-input  systems.  A 
simple  numerical  example  is  given  in  Section  9.  In  Section  10,  experimental  results  are  given 
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for  the  Deployable  Optical  Telescope  (DOT)  testbed  at  the  Air  Force  Research  Laboratory 
(AFRL).  Conclusions  are  given  in  Section  11. 

2  Problem  Formulation 

The  goal  of  frequency  domain  system  identification  is  to  obtain  a  linear  time-invariant 
model  which  approximates  the  dynamics  of  an  unknown  plant.  In  general,  the  plant  may 
have  nonlinear  time-varying  dynamics.  Plant  nonlinearities  often  induce  harmonic  responses 
to  sinusoidal  input  signals.  Consider,  for  example,  the  nonlinear  system  shown  in  Figure 
1,  which  has  asymptotically  stable  linear  dynamics  XL  with  a  static  input  nonlinearity  ATi 
and  a  static  output  nonlinearity  Af2.  The  nonlinearities  AIL  and  Af2  of  this  system  can  induce 
harmonic  responses.  In  particular,  if  AIL  and  AT2  are  analytic  functions,  then  the  nonlinearities 
can  only  induce  integer  superharmonics.  If  AR  and  AT2  are  analytic  functions,  then  they  can 
be  approximated  to  arbitrary  accuracy  by  polynomials,  that  is, 

AR(x)  —  ol 0  T  ot\x  T  •  •  •  T  oin—  \X  T  (xnx  T  ■  ■  (2.1) 

N2(.t)  =  (3q  +  / 3\x  +  •  •  •  +  Pn- \xn  1  +  f3nxn  +  •  •  •  ,  (2-2) 

where  the  real  coefficients  oto,  au, . . .  and  /?o,  /?i,  •  •  •  are  determined  by  the  Taylor  series  ex¬ 
pansions  of  Afi  and  Af2,  respectively.  Now,  we  assume  that  the  system  input  is  the  sinusoid 

u(t)  =  C\  sin  (27 rut)  +  C2  cos  (27rcat  ) .  (2.3) 

Then  the  input  to  the  linear  dynamical  system  XL  is  u&(t)  =  AR  («(£)).  Using  trigonomet¬ 
ric  identities,  u&(t)  may  be  expressed  as  a  sum  of  sinusoids  with  frequencies  a;,2u;,3u;, . . .. 
Since  XL  is  an  asymptotically  stable  linear  system  and  u&(t)  is  a  sum  of  sinusoids  with  fre¬ 
quencies  u,  2uj,  3u, . . .,  it  follows  that  the  output  yz(t)  =  XL(Af1(u(t)))  at  steady  state  is  a 
sum  of  sinusoids  with  frequencies  u,2u,3u, . . ..  Since  yn(t)  is  a  sum  of  harmonically  related 
sinusoids,  we  use  the  Taylor  series  expansion  of  N2(x)  and  trigonometric  identities  to  ex¬ 
press  y(t)  =  Af 2(yc(t))  as  a  sum  of  sinusoids  with  frequencies  o;,2n;,3cc;, . . ..  Therefore,  the 
nonlinearities  ATi  and  Af2  can  induce  integer  superharmonic  responses. 


Figure  1:  Linear  dynamic  system  with  static  nonlinearities  on  the  input  and  output. 


In  general,  a  system  may  have  nonlinear  time-varying  dynamics  that  do  not  possess  the 
block  structure  shown  in  Figure  1.  In  the  present  paper,  we  consider  plants  with  nonlinear 
time-invariant  dynamics  and  make  the  following  assumptions. 

(Al)  The  plant  has  a  globally  asymptotically  stable  equilibrium  point. 

(A2)  The  plant  is  linearizable  about  its  globally  asymptotically  stable  equilibrium. 
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(A3)  The  plant  nonlinearities  can  induce  only  integer  subharmonic  and  superharmonic  re¬ 
sponses  of  a  base  tone. 

The  objective  of  this  paper  is  to  design  an  excitation  signal  to  obtain  frequency  response 
data  for  frequency  domain  system  identification.  We  require  that  the  excitation  signal  sat¬ 
isfies  the  following  conditions. 

(i)  The  signal  must  excite  the  frequency  band  fl  =  [cnmin,  tnmax],  where  0  <  u;min  <  comax. 

(ii)  The  signal  must  be  sufficiently  rich  to  excite  frequencies  on  a  discrete  grid  with  spacing 
less  than  or  equal  to  r/  where  0  <  r]  <  cnmin. 

(in)  No  leakage  in  the  Fourier  Transform  of  the  excitation  signal  and  output  signal. 

(iv)  The  base  tone  response  of  the  plant  must  be  distinguishable  from  nonlinear  distortions 
caused  by  subharmonics  and  superharmonics. 

Without  loss  of  generality,  we  assume  that  and  are  integers. 

In  addition  to  conditions  (i)-(iv),  we  seek  an  input  signal  that  does  not  require  a  long  test 
time.  Certain  input  signals  can  require  long  test  times  to  meet  conditions  (i)-(iv).  In  fact, 
test  times  can  be  too  long  for  practical  implementation.  Thus,  the  purpose  of  this  paper  is 
to  design  an  excitation  signal  that  satisfies  (i)-(iv),  while  requiring  test  times  short  enough 
for  practical  implementation. 

In  the  following  sections,  we  examine  the  properties  of  several  common  excitation  signals. 
Specifically,  we  consider  random  noise,  burst  random  noise,  pulse-impact,  multisine,  and  sine 
dwell  excitation  signals.  For  more  information  on  these  input  signals,  see  Chapter  4  of  [8]. 
We  demonstrate  that  none  of  these  signals  satisfies  conditions  (i)-(iv),  while  being  practically 
implementable  with  respect  to  test  time. 

In  sections  3-7,  we  assume  that  the  system  is  excited  one  input  channel  at  a  time.  Thus, 
the  test  times  derived  in  sections  3-7  are  the  required  test  times  per  input  channel.  More 
time  efficient  testing  can  be  performed  by  simultaneously  exciting  all  of  the  system’s  input 
channels.  This  case  is  discussed  in  Section  8. 


3  Random  Noise,  Burst  Random  Noise,  and  Pulse-Impact 

Random  noise  excitation  requires  a  test  time  Tr  >  — —  to  achieve  the  desired  frequency 

k-’min 

resolution  and  bandwidth.  However,  test  times  Tr  »  — —  are  often  required  to  infuse  the 

k-’min 

frequency  band  with  sufficient  power,  to  obtain  a  satisfactory  signal-to-noise  ratio,  and  to 
satisfy  condition  (ii). 

Random  noise  signals  suffer  from  leakage  since  they  are  aperiodic.  However,  leakage  can 
be  alleviated  by  windowing  the  input  and  output  data  or  by  repeatedly  applying  an  input 
signal  [9].  Alternatively,  leakage  can  be  addressed  by  using  burst  random  noise  or  pulse- 
impact  excitation  signals.  Detailed  discussions  of  burst  random  noise  and  pulse-impact 
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excitation  signals  are  provided  in  [1,2].  A  burst  random  noise  signal  is 


■Ubr(l)  =  W(t)d(t), 


(3.1) 


where  w(t)  is  a  random  variable  and  d(t) 


1,  0  <  t  <  Ti 

0,  Ti  <  t  <  Tr 


A  pulse-impact  signal  is 


f  Api,  0  <  t  <  Ti 
\  0,  T\  <t  <  Tr 


(3.2) 


where  Api  e  R.  If  the  test  time  Tr  is  long  enough  so  that  the  plant’s  response  decays  to  a 
negligible  amount  by  time  Tr,  then  leakage  errors  will  be  negligible  [8].  However,  for  lightly 
damped  systems  with  low  frequency  modes,  the  plant’s  response  can  take  a  long  time  to 
decay.  In  practice,  windowing  is  often  used  with  both  burst  random  noise  and  pulse-impact 
excitation  signals  to  allow  for  shorter  test  times. 

To  make  matters  worse,  random  and  pulse  excitation  signals  do  not  provide  insight  into 
the  plant’s  nonlinearities.  Specifically,  the  plant’s  linear  response  is  indistinguishable  from 
nonlinear  distortions  since  random  noise  and  pulse  signals  simultaneously  excite  the  entire 
bandwidth  of  interest.  Thus  random  noise,  burst  random  noise,  and  pulse-impact  signals  do 
not  satisfy  condition  (tv). 


4  Multisine 

A  multisine  signal  is  a  summation  of  sinusoids  with  varying  frequencies  and  phases.  A 
Schroeder-phased  multisine  is  the  sum  of  sinusoids  where  the  phases  are  selected  to  obtain 

low  peak-to-peak  amplitudes  [3].  Define  N  =  ^max~^min  +1  and  consider  the  Schroeder-phased 
multisine 


N 

=  E  Am  cos  (2iTUkt  +  4> k) , 

k=  i 


where  Am  e  R  and,  for  k  —  1, . . . ,  N, 

A  —k(k  —  l)n 
<Pk  =  - jy - . 

=  ^min  T  T]  {k  1)  . 


(4.1) 


(4.2) 

(4.3) 


The  Schroeder  phases  (pi, . . . ,  (j>N  have  been  used  for  simplicity.  However,  the  phases  may  be 
selected  to  obtain  other  desired  properties  of  the  input  signal.  For  example,  the  phases  may 
be  selected  to  minimize  the  crest-factor  [10].  Note  that  the  multisine  (4.1)  has  a  flat  power 
spectrum  over  Q. 

The  multisine  is  periodic  with  period  A  Therefore,  leakage  is  eliminated  if 

Tm  =  lm  f  A tm,  (4-4) 

V 
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where  lm  E  N  is  the  number  of  periods  used  for  averaging  and  A tm  >  0  is  the  time  required 
for  the  system’s  response  to  um(t)  to  converge  to  steady  state.  Multisines  can  be  efficient 
with  respect  to  test  time  since  they  simultaneously  excite  the  entire  frequency  band  0  and 
require  only  one  settling  period  A tm. 

However,  simultaneous  excitation  of  the  entire  bandwidth  0  is  undesirable  for  dealing 
with  nonlinearities.  It  is  impossible  to  separate  the  linear  response  from  harmonic  responses 
caused  by  plant  nonlinearities,  and  thus  multisine  excitation  signals  do  not  satisfy  condition 
(iv). 


5  Sine  Dwell 

A  sine  dwell  excitation  signal  is  a  sequence  of  sinusoidal  excitation  signals,  each  with 
different  frequency.  Define  N  =  CJmax"a-’min  4- 1  and,  for  k  =  1 , ,N,  consider  the  sinusoids 

us,k(t)  =  AStk  sin  (2 nujkt) ,  (5.1) 

where  for  k  —  1, . . . ,  N,  As,k  e  M  and  u>k  =  u;min  +  77  ( k  —  1).  Sine  dwell  testing  is  performed 
by  conducting  the  sequence  of  N  sinusoidal  inputs  uS:i(t), . . .  ,  ws,iv(£)-  This  populates  the 
frequency  band  D  with  frequency  spacing  r/. 

The  sinusoids  uS:i(t), . . .  ,uSiAr(t)  are  periodic  with  periods  T-, . .  ^  A-,  respectively.  For 
k  =  1, . . . ,  N,  let  TS:k  be  the  measurement  time  associated  with  uStk(t).  Leakage  is  eliminated 
by  letting  TS)k  =  ls,k^  where>  for  k  =  1, . . . ,  N,  ls  k  G  N  is  the  number  of  periods  used  for 
averaging. 

LInlikc  the  other  excitation  signals  we  have  discussed  thus  far,  a  sine  dwell  excitation 
signal  allows  the  harmonic  response  to  be  separated  from  the  base  tone  response.  For 
k  =  1 , ,N,  the  sinusoid  us,k(t)  is  used  to  estimate  the  base  tone  response  at  frequency 
uikl  but  the  harmonic  response  can  be  distinguished  from  the  base  tone  response  since  the 
harmonic  responses  due  to  nonlinearities  occur  at  frequencies  other  than  u>k.  In  fact,  sine 
dwell  testing  provides  detailed  information  of  nonlinear  distortions. 

However,  the  major  drawback  to  sine  dwell  testing  is  lengthy  test  time.  The  test  time 
required  for  sine  dwell  is 


N 


N 


Ts  —  Ts,k  +  A tS)k  —  y]  ZSjfc  b  AfSjfc, 


k=  1 


fc=l 


Uk 


(5.2) 


where,  for  k  =  1, . . . ,  N,  A tSjk  >  0  is  the  time  required  for  the  system’s  response  to  the 
sinusoid  uSjk(t)  to  settle.  Sine  dwell  testing  can  require  significantly  longer  test  times  than 
multisine  testing. 

Table  1  summarizes  the  properties  of  the  excitation  signals  discussed  in  sections  3-5.  Only 
sine  dwell  allows  the  linear  frequency  response  to  be  separated  from  nonlinear  distortions. 
However,  sine  dwell  can  require  lengthy  test  times. 
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Excites  Q 

Frequency 
Spacing  r] 

No  Leakage 

Removes  Distortions 
due  to  Harmonics 

Test  Time  T 

Random  Noise 

/ 

Often  requires 

T»  zr— 

k-’mm 

Requires  windowing 

No 

Condition  (ii) 
necessitates  large  T 

Burst  Random  Noise 

/ 

Often  requires 

T»  zr~ 

Requires  T  >>  — - — 

wmin 

and/or  windowing 

No 

Conditions  (ii)-(iii) 
necessitates  large  T 

Pulse-Impact 

/ 

Often  requires 

T»  zr~ 

^mm 

Requires  T  >>  — - — 

wmin 

and/or  windowing 

No 

Conditions  (ii)-(iii) 
necessitates  large  T 

Multisine 

S 

s 

No 

More  time  efficient 

Sine  Dwell 

s 

s 

V 

Can  require  large  T 

Table  1:  Properties  of  random  noise,  burst  random  noise,  pulse-impact,  multisine,  and  sine 
dwell  excitation  signals. 


6  Test  Times  for  Multisine  and  Sine  Dwell  Excitation 


In  this  section,  we  compare  the  test  times  required  for  multisine  and  sine  dwell  excitation 
signals.  Specifically,  we  compare  the  test  times  required  to  obtain  frequency  response  data 
with  equivalent  signal-to-noise  ratios.  The  quality  of  system  identification  data  is  most  accu¬ 
rately  described  by  the  signal-to-noise  ratios  of  the  input  and  output  signals.  However,  input 
signals  are  usually  assumed  to  be  noise-free,  and  the  output  signal-to-noise  ratio  depends  on 
the  output  noise,  the  input  power  spectral  density,  and  the  system  dynamics.  If  the  system 
is  linear,  then  the  output  power  spectral  density  is  the  product  of  the  input  power  spectral 
density  and  the  magnitude  of  the  system’s  transfer  function.  Although  this  paper  considers 
nonlinear  systems,  if  input  amplitudes  are  small,  then  the  output  power  spectral  density 
is  approximated  by  the  product  of  input  power  spectral  density  and  the  magnitude  of  the 
linearized  transfer  function.  We  use  the  signal-to-noise  ratio  obtained  from  the  approximate 
output  power  spectral  density  as  a  measure  of  data  quality. 


Let  Fs  denote  the  sampling  frequency  of  the  multisine  and  sine  dwell  excitation  signals. 
The  discrete  Fourier  transform  of  one  period  of  the  multisine  um(t)  is 


Um(k)  = 


(6.1) 


and  the  discrete  Fourier  transform  of  an  excitation  signal  consisting  one  period  of  each  sine 
dwell  . . .  ,uStisr(t)  is 


Ua(k) 


[FsA^k 

V  ^ k  2 


(6.2) 


Assumption  (Al)  implies  that  the  excitation  signals  are  for  linearizable  systems,  so  let 
G(z )  be  the  discrete  transfer  function  matrix  of  the  linearized  and  sampled  system.  Then 
the  approximate  output  power  spectral  densities  obtained  from  multisine  and  sine  dwell 
excitation  are 


Ym(k)  = 

Ys(k)  = 


\\G(e™)\\ 


l|G(e^)||, 

IIGV^II, 


(6.3) 

(6.4) 
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y/N 


Figure  2:  The  peak  signal  value  um  of  a  Schroeder-phase  multisine  versus  yj N.  Various 
values  of  77  are  plotted  where  u;min  =  rj. 


respectively. 


Input  amplitudes  for  the  mutlisine  and  sine  dwell  excitation  signals  must  be  normalized 
to  compare  testing  times,  since  linear  approximations  are  only  valid  for  small  amplitudes 
and  actuator  saturation  imposes  practical  limitations.  If  we  design  multisine  and  sine  dwell 
excitation  signals  such  that  Am  =  A, ;1  =  •  •  •  =  A,jv  =  A  where  A  >  0,  then  the  excitation 
signals  will  not  have  equal  peak  values.  In  fact,  the  peak  value  of  a  multisine  increases  with 
the  number  of  sinusoids  N,  see  Figure  2. 

Define  the  peak  value  um  =  max0<t<i 

multisine  (4.1).  Figure  2  provides  plots  of  um  versus  yN  for  various  values  of  rj  where  u;min  = 
77.  The  numerical  results  indicate  that  the  peak  value  of  the  Schroeder-phased  multisine  (4.1) 
is  proportional  to  y/N.  Furthermore,  the  constant  of  proportionality  is  independent  of  the 
frequency  spacing  77. 


Ylk=i  cos  {ZnuJkt  +  A) 


of  the  Schroeder-phased 


Let  A.i  =  •  •  •  =  A,jv  =  A  and  Am  =  ^=,  where  A  >  0  is  the  desired  peak  signal 
value.  We  use  the  normalization  Am  =  which  is  convenient  for  the  subsequent  analysis, 

even  though  the  constant  of  proportionality  between  um  and  V N  is  not  exactly  one.  The 
approximate  output  power  spectral  densities  satisfy 


(6.5) 

(6.6) 


Let  Ny(k)  denote  the  discrete  Fourier  transform  of  the  output  noise,  then  the  output 
signal-to- noise  ratio  for  a  single  period  of  the  multisine  excitation  um(t)  is 


Ym(k) 

Ny(k) 


1  A^F~s 
y/rjN  2Ny(k) 


HG(e^)||, 


(6.7) 
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and  the  output  signal-to-noise  ratio  for  an  excitation  signal  with  a  single  period  of  each  sine 
dwell  uS:i(t), . . .  ,uS:isr(t )  is 


YB(k) 

Ny(k) 


1  A^FS 
y/u^2Ny{k) 


||G(e^)||. 


(6.8) 


The  test  times  required  for  the  multisine  and  sine  dwell  excitations  are  given  by  (4.4) 
and  (5.2),  respectively.  However,  the  number  of  multisine  periods  lm  and  sine  dwell  periods 
4,i,  ■  •  • ,  4,jv  used  for  averaging  should  be  chosen  to  ensure  that  the  output  signals  have 
comparable  signal-to-noise  ratios. 

We  assume  that  Ny(k)  is  filtered  white  noise,  independent  of  both  Ym(k )  and  Ys(k),  with 
a  finite  second  moment.  Under  these  mild  assumptions  averaging  the  output  data  over  l  G  N 
periods  reduces  the  noise  proportional  to  ^  [8].  Let  SNRm(/c)  denote  the  output  signal-to- 
noise  ratio  of  the  multisine  excitation  data  when  lm  averages  are  performed,  and  let  SNRs(A;) 
denote  the  output  signal-to-noise  ratio  of  the  sine  dwell  excitation  data  when  4,i,  •  •  •  ,4 ,n 
averages  are  performed.  Therefore,  the  signal-to-noise  ratios  satisfy 

SNMfc)  .X  Eu  llG(e*‘l‘)ll  ’  <6'9> 

and 

SNR.WOC  (6.10) 


If  4,fc  =  /yy  and  lm  =  IN,  then  it  follows  from  (6.9)  and  (6.10)  that  the  output  signals 
obtained  from  the  multisine  and  sine  dwell  excitations  have  comparable  signal-to-noise  ratios. 
To  obtain  these  comparable  signal-to-noise  levels,  the  multisine  excitation  signal  requires 


Tm  —  4n — b  Atn 
V 


N — b  A tm, 
V 


and  sine  dwell  excitation  requires 


N  i  iV  /  \  i  l  N 

Ts  =  4,fc - b  A ts  k  =  ( l —  ) - b  A ts  k  =  N — b  A ts  k- 


(6.11) 


(6.12) 


It  follows  from  (6.11)  and  (6.12)  that  sine  dwell  testing  can  take  significantly  longer  than 
multisine  testing.  Specifically,  sine  dwell  testing  requires  additional  time  for  the  system  to 
settle  each  time  the  sine  dwell  frequency  is  changed.  For  linear  systems,  practical  settling 
times  A tm  and  AfS;1, . . . ,  Ats,jv  can  be  estimated  from  the  system’s  time  constant.  Further¬ 
more,  the  time  constant  of  a  linear  system  is  independent  of  input  signal.  Therefore,  in 
practice,  it  is  often  assumed  that  settling  time  is  independent  of  the  harmonic  content  of  the 
excitation  signal. 

If  we  assume  that  the  settling  times  A fm  and  AtSti, . . .  ,Ats,N  are  independent  of  exci¬ 
tation  signal,  then  sine  dwell  testing  requires  ( N  —  1)  A tm  time  units  longer  than  multisine 


testing.  If.  in  addition,  the  frequency  band  of  interest  is  large  and/or  the  frequency  spacing 
is  small,  then  N  will  be  large  and  sine  dwell  testing  will  take  significantly  longer  than  mul¬ 
tisine  testing.  Furthermore,  if  the  system  is  lightly  damped,  then  settling  times  can  be  very 
long,  exacerbating  the  disparity  in  test  times. 

For  example,  consider  a  single-input  system  and  design  input  signals  to  excite  frequencies 
from  um in  =  0.1  Hz  to  o;max  =  100  Hz  with  a  frequency  spacing  of  rj  =  0.1  Hz.  Assume  that 
the  system  requires  5  seconds  to  settle.  The  multisine  signal  (4.1)  requires 

Tm  =  N — b  A tm  =  10,  000/  +  5  sec,  (6.13) 

V 

and  the  sine  dwell  signal  (5.1)  requires 

/  N 

Tsk  —  N-  +  W  A tsk  =  10,  000/  +  5,  000  sec.  (6.14) 

77  ^ 

1  k= 1 

Thus,  sine  dwell  testing  requires  an  additional  4,995  seconds  (or  1  hour  23  minutes  15 
seconds).  If  /  is  small,  then  4,995  seconds  constitutes  a  significant  portion  of  the  total 
testing  time.  Furthermore,  if  the  frequency  band  fl  is  increased  or  the  frequency  spacing  r) 
is  decreased,  then  the  disparity  between  the  sine  dwell  test  time  and  the  multisine  test  time 
will  increase. 


7  Sequential  Multisines 

In  this  section,  we  present  a  novel  excitation  signal  consisting  of  a  sequence  of  multi¬ 
sines  that  are  designed  to  retain  the  benefits  of  sine  dwell  signals  while  being  more  time 
efficient.  Specifically,  we  design  a  sequence  of  multisines  that  allow  the  linear  response  to  be 
distinguished  from  nonlinear  distortions. 

Let  M  e  N  be  the  least  positive  integer  such  that  M  >  log2  In  general,  M  <  N 

^min 

and  if  the  bandwidth  of  interest  is  large  and/or  r]  is  small,  then  M  «  N . 

Now,  we  divide  the  frequency  band  H  into  M  segments.  For  i  =  1, . . . ,  M,  define  the 
frequency  bands 

a  [-A  •  ■  i )  i  (7.1 ) 

where,  for  7  =  1,...,  M, 

=  2‘-1aW,  (7.2) 

and  Qm+i  =  +  ^  Note  that,  for  i  =  1, . . . ,  M,  does  not  contain  integer  subharmonics 

or  superharmonics  of  the  base  frequencies  contained  in  H*. 

Now,  we  consider  a  sequence  of  M  Schroeder-phased  multisines.  For  i  —  1, ,  M,  define 
Ni  =  u)i+1~u)i  and  consider 

L  ^ 

Ni 

^sm  ,i  «  =  £  -^sm  ,i  COS  (27 TLUijt  +  4>ij)  ,  (7.3) 

j= 1 
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(7.4) 

(7.5) 


where,  for  i  —  1, ... ,  M,  Asnvj  G  M  and,  for  j  —  1, . . . ,  Nl} 


A  -j(;j  -  l)?r 


'h3 


Nj 


A 


Ujid  =uJi  +  r)  ( j  -  1) 


If,  for  example,  cjmin  =  rj,  then  Nx  —  1,  N2  =  2,  N3  —  4,.  jV4  —  8, . . . ,  NM_X  =  2M_2  and 
A4v/  =  2m~2  +  +  1.  Again,  the  Schroeder  phases  have  been  chosen  for  simplicity,  but 

the  phases  may  be  chosen  to  satisfy  other  excitation  signal  criteria. 

For  i  —  1, ,  M,  the  multisine  usm_i(t)  excites  frequencies  in  the  band  but  does  not 
excite  frequencies  outside  of  the  band  44 .  For  i  —  1, ... ,  M,  let  Tsmii  be  the  measurement  time 
associated  with  usm^{t).  For  i  —  1, . . . ,  M,  letting  Tsul)i  =  4m, where  4 m>i  G  N  eliminates 
leakage.  Then  the  test  time  required  for  the  sequential  multisine  signal  is  given  by 


M 

Tsm  =  ^2  ,i 
2—1 


1  M 

Atsm  i  /  4m A 

7] 

1  7=1 


+  At 


sm,D 


(7.6) 


where,  for  i  —  1, . . . ,  M,  Aism,i  >  0  is  the  time  required  for  the  system’s  response  to  usmj(t) 
to  settle  to  harmonic  steady  state. 

Now,  we  compare  the  sequential  multisine  test  time  with  the  multisine  and  sine  dwell 
test  times.  For  i  =  1 , . . . ,  M,  let  Asulti  =  -^=,  thus  normalizing  the  peak  signal  values 
to  approximately  A  >  0.  The  number  of  sequential  multisine  periods  4m, l,  ■  •  • ,  4m, m  used 
for  averaging  must  be  chosen  to  ensure  that  the  output  signal  obtained  from  the  sequential 
multisine  excitation  has  a  signal-to-noise  ratio  comparable  to  those  obtain  from  the  multisine 
and  sine  dwell  excitations. 


For  i  =  1, ... ,  M ,  the  discrete  Fourier  transform  of  a  signal  consisting  of  one  period  of 
each  sequential  multisine  usm,i, . . . ,  usm,M  is 


44m  (4) 


[Fs  A 

V  V  2  V3VT 


\/k  G  [A4_i,  Ni] , 


(7.7) 


where  N0  =  1.  Then,  for  i  —  1, . . . ,  M,  the  approximate  output  power  spectral  density  is 

y„(fc)i-A=A5||G(e^.)||,  VA:  6  [«.,,«],  (7.8) 

which  implies  the  output  signal-to-noise  ratio  is 


Ysm(k)  _  1  A^/F~s 

A Ty(k)  \/r)Ni  2Ny(k) 


||G(eJWfc)||l  V4  G  [A4j_1;  A4] 


(7.9) 


Let  SNRsm(4)  denote  the  output  signal-to-noise  ratio  of  the  sequential  multisine  excitation 
data  when  for  *  =  1, ... ,  M,  4m, i  averages  are  performed.  Then,  for  i  —  1, . . . ,  M, 


SNR,sm(£;)  oc  s/Ts 


AVZ 


\\G(e>^)\\,  V4  G  [A4.-i,  Ni\  ■ 


(7.10) 
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If,  for  i  =  1, . . . ,  M,  /sm =  /A*,  then  SNRsm(4)  is  comparable  to  SNRm(/c)  and  SNRs(4).  To 
obtain  this  comparable  signal-to-noise  level,  the  sequential  multisine  excitation  requires 


M  i  M  i  M 

Tsm  =  Tsm  j  +  A ZAj  +  A4m,j  =  A- - h  A4m,i-  (7-11) 


2=1 


^  AT 


v 


2=1 


It  follows  from  (6.11),  (6.12),  and  (7.11)  that  sequential  multisine  testing  is  less  time 
efficient  than  single  multisine  testing,  but  sequential  multisine  testing  is  more  time  efficient 
than  sine  dwell  testing.  In  fact,  if  the  system’s  settling  times  A tm,  Ats,i, . . . ,  A4,jv,  and 
A4,i,  •  •  • ,  A4,jv  are  equal,  then  sine  dwell  testing  takes  (A  —  M)Atm  more  time  units  than 
sequential  multisine  testing.  Furthermore,  A  increases  quicker  than  M  as  the  number  of 
discrete  frequency  points  increases.  Specifically,  A  increases  linearly  with  the  number  of 
discrete  frequencies  whereas  M  is  proportional  to  the  logarithm  of  the  number  of  discrete 
frequencies.  Therefore,  A  —  M  may  be  very  large. 

Consider  again  the  example  given  in  Section  6.  The  objective  is  to  design  an  input  signal 
that  excites  frequencies  from  =  0.1  Hz  to  cnmax  =  100  Hz  with  a  discrete-frequency 
spacing  of  r/  =  0.1  Hz.  We  assume  that  for  all  harmonic  excitations,  the  system  requires  5 
seconds  to  settle.  Recall  that  the  a  single  multisine  requires  a  test  time  of  Tm  =  10,  000/  +  5 
sec,  and  a  sine  dwell  input  requires  a  test  time  of  Ts  =  10,  000/  +  5, 000  sec.  In  contrast,  the 
sequential  multisine  (7.3)  requires  a  test  time 


/ 


M 


Tsm  =  A — b  A4m  i  =  10,  000/  +  50  sec. 

77 

1  i= 1 


(7.12) 


Thus  sine  dwell  excitation  requires  4,950  seconds  (or  1  hour  22  minutes  30  seconds)  longer 
than  sequential  multisine  excitation. 

For  7  =  1,...,  M,  usm  i  (/)  is  used  to  estimate  the  base  tone  response  over  the  frequency 
band  f 4,  but  the  harmonic  response  due  to  nonlinearities  occur  at  frequencies  not  contained 
in  the  fi?.  Thus,  the  linear  response  over  0,:  is  distinguishable  from  the  nonlinear  distortions 
that  occur  at  frequencies  outside  of  Qt.  Thus,  the  sequential  multisine  satisfies  conditions 
(i)-(iv)  while  being  less  time  consuming  than  sine  dwell  testing. 


8  Sequential  Multisine  Excitation  for  Multi-Input  Systems 

In  the  previous  sections,  we  designed  multisine,  sine  dwell,  and  sequential  multisine 
excitation  signals  assuming  that  the  system  is  excited  one  input  at  a  time,  that  is,  data  is 
first  taken  by  exciting  input  channel  1,  then  by  exciting  input  channel  2,  and  continuing  until 
the  system  has  been  excited  through  all  input  channels.  Now,  we  assume  that  the  system 
has  p  inputs.  Then  the  test  times  given  by  (6.11),  (6.12),  and  (7.11)  are  multiplied  by  p  to 
yield  the  total  test  time. 

For  a  sine  dwell  excitation  signal,  it  is  clear  that  exciting  the  system  using  only  one 
input  channel  at  a  time  is  not  time  efficient.  In  fact,  sine  dwell  testing  can  be  performed 
by  simultaneously  exciting  all  input  channels.  To  ensure  that  we  can  distinguish  between 
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the  different  input-to-output  frequency  responses,  we  require  that  no  two  input  channels  are 
simultaneously  excited  at  the  same  frequency.  Similarly,  to  ensure  that  the  base  tone  response 
is  distinguishable  from  the  nonlinear  distortions,  we  require  that  no  two  input  channels  be 
simultaneously  excited  at  any  harmonic  of  the  same  base  frequency.  Additionally,  data  must 
be  taken  for  an  integer  number  of  periods  at  each  frequency  to  avoid  leakage  in  the  Fourier 
Transform.  Following  these  rules,  the  test  time  for  sine  dwell  excitation  is  independent  of 
the  number  of  inputs. 

Similarly,  sequential  multisine  excitation  signals  can  be  designed  so  that  all  input  channels 
may  be  excited  simultaneously,  yielding  a  test  time  that  is  independent  of  p.  We  consider  the 

frequency  bands  f21; . . .  ,Om  given  by  (7.3),  and  defined  the  finer  frequency  spacing  fj  =  A 
For  each  input  channel,  we  consider  a  sequence  of  M  Schroeder-phased  multisines.  For 
i  —  1, . . . ,  M,  consider 

Ni 

Ani.y.  (4)  ^  ^  ASm,i  COS  (2tT  ( ^i,j )  t  T  (pi  J  )  , 

3= 1 
Nt 

U2sm ,i(f)  =  X  Asm’{  C°S  (2n  +  V)  *  +  , 

3= 1 
Ni 

«L ,i(t)  =  X  cos  (2?r  (W*J  +  2fi ) 1  +  M  > 

3= 1 

:  (8.4) 

Nt 

<n,;  W  =  cos  (2?r  +  (P  “  x)  V)  t  +  <pi,j)  ■  (8.5) 

j=i 


(8.1) 

(8.2) 

(8.3) 


In  this  notation,  input  channel  1  is  excited  with  the  sequence  of  multisines  u)m  i(t). . . . ,  ulmM(t), 
input  channel  2  is  excited  with  the  sequence  of  multisines  v‘lm  l(t), . . .  and  input 

channel  p  is  excited  with  the  sequence  of  multisines  v!P,m  ,  (f), . . .  ,ups 

For  i  =  1, ,  M,  the  multisines  ulm  i(t), . . . ,  are  periodic  with  period  L  There¬ 

fore,  leakage  is  eliminated  by  letting  the  test  time  associated  with  each  multisine  in  the 
sequence  be  Tsrn  t  =  lsm/(P  where  lsm,t  €  N.  Since  all  input  channels  are  simultaneously 
excited,  the  required  test  time  is  given  by 


M 


T  =  - 

■J-  sm 


V 

1  i=  1 


(8.6) 


where,  for  i  —  1, . . . ,  M,  A tsm,j  is  the  time  required  for  the  system’s  response  to  ulm  i(t), . . . , 
to  settle. 

If  data  is  collected  for  at  least  4  seconds  and  the  settling  time  is  independent  of  the 
excitation  signal,  then  the  test  time  required  for  sequential  mutlisine  excitation  does  not 
depend  on  the  number  of  inputs.  That  is,  the  system  can  be  excited  though  all  channels 
simultaneously. 
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y{t)  u(t) 


Figure  3:  Two-degree-of-freedom  spring-mass-damper  with  nonlinear  stiffness. 

9  Two-Degree-Of-Freedom  Spring-Mass-Damper  with  Cubic  Stiff¬ 
ness 


Consider  the  two-degree-of-freedom  spring-mass-damper  with  nonlinear  stiffness  shown 
in  Figure  3.  The  dynamics  of  the  system  are 


Mq  +  Cq  +  K{q)  =  bu, 


where 


M 

C 

K(q) 

q 


A 

rrii 

h 

A 

'  o ' 

m2 

i  u 

i 

A 

Cl  +  c2  - 

-C2 

c2 

c2 

) 

A 

ki(q i)  -  k2(q2 

-<?i)  ’ 

l  k2(q 2 

-Qi) 

j 

1T 

A  | 

Qi  Q2 

(9.1) 

(9.2) 

(9.3) 

(9.4) 

(9.5) 


Assume  that  the  nonlinear  stiffness  ki(-)  and  A?2 (•)  are  given  by  the  cubic  stiffness  terms 

k\{x)  =  k\x  +  k\x3  and  k2{x)  =  k2x  +  k2x3.  The  masses  are  m\  =  m2  =  0.1  kg.  The 
damping  coefficients  are  C\  =  0.1  kg/sec  and  ci  =  0.05  kg/sec.  The  spring  constants  are 
k\  =  56  kg/sec2,  k2  =  90  kg/sec2,  and  k\  —  k2  —  50  kg/sec2/m2.  The  system  (9.1)-(9.5)  is 
linearizable  about  the  origin,  which  is  an  asymptotically  stable  equilibrium.  Let  the  output 
y  be  the  position  of  the  first  mass,  that  is  y  —  q\ .  The  linearized  system  from  input  u  to 
output  y  is 


u  =  G(s)y, 


(9.6) 


where 


5s  +  9000 

s4  +  2s3  +  2360.5s2  +  1180s  +  504000 
5s  +  9000 

(s2  +  0.374s  +  237.464) (s2  +  1.626s  +  2122.4) ' 


(9.7) 
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The  system  has  modes  at  approximately  2.45  Hz  and  7.33  Hz.  Now  let  us  consider  excitation 
signals  designed  for  frequency  domain  system  identification  over  the  frequency  band  cnmin  = 
0.1  Hz  and  cnmax  =  25  Hz  with  a  frequency  spacing  of  rj  =  0.1  Hz.  The  peak  value  of  the 
excitation  signal  must  be  normalized  to  approximately  A  =  0.25. 

We  consider  the  following  three  excitation  signals. 

(1)  The  single  multisine  given  by  (4.1)  where  N  =  250  and  Am  =  N. 

(2)  The  sine  dwell  signal  (5.1)  where  N  =  250  and  HS)i, . . . ,  ASt N  =  0.25  N. 

(3)  The  sequential  multisine  (7.3)  where  M  —  8  >  log2(^ai2:2L)  and,  for  i  =  1 , . . . ,  M, 

Asm,i  =  ^  where  Nt  =  Wi+1~Wi ,  Therefore,  Hsma  =  0.25  N,  Hsni)2  «  0.1768  N, 
nsm,3  =  0.125  N,  Hsm,4  ~  0.0884  N,  Hsm)5  =  0.0625  N,  Asmfi  «  0.0442  N,  Hsm>7  «  0.0313 
N,  and  Hsmi8  ~  0.0221  N. 


To  obtain  equivalent  signal-to-noise  ratios,  the  single  multisine  excitation  is  averaged 
over  lm  =  N  periods,  for  k  —  1, . . . ,  N,  the  sine  dwell  excitations  are  averaged  over  ls  k  =  ^ 
periods,  and,  for  i  —  1, . . . ,  M,  the  sequential  multisine  excitations  are  averaged  over  4m, i  = 
Nt  periods. 


To  illustrate  the  effects  of  nonlinearity  induced  harmonics,  we  consider  the  case  where 
there  is  no  noise  on  the  input  or  the  output.  We  excite  the  system  (9.1)-(9.5)  using  the  single 
multisine  excitation.  The  system  is  allowed  A tm  =  10  sec  to  settle.  Then  data  is  collected 
for  lm  =  N  =  250  periods  of  the  single  multisine,  yielding 

Tm  =  lm — b  A tm  =  2, 510  sec  =  41  min  50  sec.  (9.8) 

V 

The  transfer  function  (9.7)  of  the  linearized  system  and  the  frequency  response  data  obtain 
using  the  single  multisine  are  shown  in  Figure  4.  Input-output  data  is  collected  for  2,500  sec 
at  a  sample  rate  of  500  Hz.  To  obtain  the  frequency  response  data,  we  use  Welch’s  averaged 
periodogram  method,  averaging  250  square  windows  of  5,000  data  points  to  obtain  5,000 
fast  Fourier  transform  points. 


Next,  we  excite  the  system  (9.1)-(9.5)  with  the  sine  dwell  signal.  The  system  is  allowed 
AfSji,  •  •  • ,  A 4,tv  =  10  sec  to  settle  at  each  frequency.  For  k  —  1, ... . . ,  N,  we  collect  data  for 
4,fc  =  qy  =  4  periods  of  the  sinusoid  u8tk(t).  Thus, 

N  1 

Ts  =  4  k - b  A ts  k  =  5,  000  sec  =  1  hr  23  min  20  sec.  (9.9) 

The  transfer  function  (9.7)  of  the  linearized  system  and  the  frequency  response  data  obtain 
using  the  sine  dwell  input  are  shown  in  Figure  5. 


Finally,  the  system  (9.1)-(9.5)  is  excited  with  the  sequential  multisine,  and  allowed 
AtmS)i, . . . ,  AtmSjM  =  10  sec  to  settle.  For  i  =  1 , . . . ,  M,  we  collect  data  for  4m, *  = 
periods  of  'usm,i(f),  yielding 


M 


Tsm  =  4m  i — b  A 4m  i  =  2,580  sec  =  43  min. 

^ J  ’  i] 


i— 1 


(9.10) 
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Figure  4:  The  transfer  function  of  the  linearized  system  (solid)  and  the  frequency  response 
data  obtain  using  the  single  multisine  excitation  (dashed). 


Figure  5:  The  transfer  function  of  the  linearized  system  (solid)  and  the  frequency  response 
data  obtain  using  the  sine  dwell  excitation  (dashed). 
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Figure  6:  The  transfer  function  of  the  linearized  system  (solid)  and  the  frequency  response 
data  obtain  using  the  sequential  multisine  excitation  (dashed). 


The  transfer  function  (9.7)  of  the  linearized  system  and  the  frequency  response  data  obtain 
using  the  sequential  multisine  are  shown  in  Figure  6. 

Figures  4-6  illustrate  that  the  frequency  response  data  obtained  from  the  single  multisine 
excitation  provides  the  poorest  approximation  the  linearized  dynamics.  Nonlinear  distortions 
in  the  frequency  response  occur  above  the  7.33  Hz  mode  and  are  a  result  of  the  third  harmonic 
of  the  2.45  Hz  mode.  The  frequency  response  data  obtained  from  the  sine  dwell  excitation, 
shown  in  Figure  5,  and  the  frequency  response  data  obtained  from  the  sequential  multisine 
excitation,  shown  in  Figure  6,  both  provide  good  approximations  of  the  linearized  dynamics. 
Furthermore,  if  the  input  and  output  data  obtained  from  the  sine  dwell  and  sequential 
multisine  excitations  were  corrupted  by  the  same  noise  sequence,  then  the  associated  noise 
levels  of  the  two  sets  of  frequency  response  data  would  be  comparable.  However,  the  sine 
dwell  testing  takes  40  minutes  20  seconds  longer  than  the  sequential  multisine  testing. 


10  Deployable  Optical  Telescope 

The  Deployable  Optical  Telescope  (DOT)  is  a  space-traceable  sparse-aperture  telescope 
developed  by  the  Air  Force  Research  Laboratory,  see  Figure  7.  The  weight  and  volume 
restrictions  imposed  by  launch  vehicles  are  the  main  impediments  to  fielding  large  aperture 
space  telescopes.  To  alleviate  volume  restrictions,  DOT’s  primary  mirror  segments  and 
secondary  mirror  tower  are  deployable.  Specifically,  the  primary  mirror  segments  rotate  into 
place  via  hinges  and  lock  with  latches.  The  secondary  mirror  tower  slides  into  place  and 
locks  with  latches. 

A  precise  alignment  of  DOT’s  optical  path  is  required  to  meet  imaging  requirements. 
Therefore,  DOT  requires  active  feedback  control  to  reject  image  path  disturbances  caused 
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Figure  7:  The  Deployable  Optical  Telescope  (DOT)  testbed. 


by,  for  example,  thermal  variations  and  the  motion  of  the  spacecraft’s  reaction  wheels.  Each 
primary  mirror  segment  is  actuated  with  three  degrees  of  freedom,  namely  tip,  tilt,  and 
piston.  The  secondary  tower  is  controlled  using  a  piezo  patch  bonded  to  the  tower.  DOT  uses 
laser  metrology  to  provide  optical-precision  measurements  for  feedback  control.  Specifically, 
the  laser  metrology  system  provides  tip,  tilt,  and  piston  measurements  for  each  of  the  three 
primary  mirror  segments.  In  sum,  DOT  has  10  control  inputs  and  9  feedback  measurements. 
For  more  information  on  DOT  see  [11-13]. 

High-fidelity  linear  dynamic  models  of  DOT  are  required  to  design  active  feedback  con¬ 
trollers  that  reject  disturbances  to  DOT’s  optical  path.  However,  obtaining  input-output 
data  for  system  identification  of  DOT  is  difficult  because  system  nonlinearities  induce  har¬ 
monic  responses  [14] .  Sine  dwell  excitation  signals  could  be  used  to  obtain  frequency  response 
data  of  DOT’s  base  tone  responses,  but  this  would  be  very  time  consuming  because  perfor¬ 
mance  objectives  require  frequency  response  data  over  a  large  frequency  bandwidth  (0  to 
1000  Hz),  and  DOT  has  densely  spaced  modes,  requiring  frequency  response  data  with  a 
small  frequency  spacing.  Therefore,  system  identification  data  may  be  obtained  using  the 
sequential  multisine  excitation  signal. 

Now,  we  compare  the  frequency  response  data  obtained  from  DOT  when  using  single 
multisine,  sine  dwell,  and  sequential  multisine  excitation  signals.  DOT  is  excited  through 
the  first  actuator  on  the  first  primary  mirror  segment  and  the  piston  position  of  the  first 
primary  mirror  segment  is  measured.  Thus,  we  obtain  frequency  response  data  from  the 
first  actuator  of  the  first  mirror  to  the  piston  position  of  the  first  mirror.  To  compare  the 
exictation  signals,  frequency  response  data  is  obtained  over  the  frequency  band  ojm\ n  =  0.5 
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Hz  and  o;max  =  90  Hz  with  a  frequency  spacing  of  rj  =  0.5  Hz.  The  peak  value  of  the 
excitation  signal  must  be  normalized  to  approximately  A  =  0.02.  We  consider  the  following 
three  excitation  signals. 


(1)  The  single  multisine  given  by  (4.1)  where  N  =  250  and  Am  =  =  A2L  N. 

(2)  The  sine  dwell  signal  (5.1)  where  N  =  180  and  HS)1, . . . ,  ASyN  =  0.02  N. 

(3)  The  sequential  multisine  (7.3)  where  M  =  8  >  log2(^^SMi2L)  and,  for  i  = 
Am,  =  7=  where  Nt  = 


To  obtain  equivalent  signal-to-noise  ratios,  the  single  multisine  excitation  is  averaged 
over  lm  =  360  periods,  for  k  =  the  sine  dwell  excitations  are  averaged  over  4,fc  = 

2^  periods,  and,  for  %  =  1,...,M,  the  sequential  multisine  excitations  are  averaged  over 
4m,*  =  2 Ni  periods. 

For  all  excitation  signals,  the  system  is  allowed  A tm  =  Ats,i, . . . ,  A =  AimS)1, . . . ,  A tmS)M 
10  sec  to  settle.  Therefore,  the  single  multisine,  sine  dwell,  and  sequential  multisine  test  times 
are 


T 

1  m 

T 

-*■  s 


lm — h  A tm  =  730  sec  =  12  min  10  sec, 
V 


N 


1 


E4fc - h  A  ts  k  =  2,  520  sec  =  42  min, 


k= 1 
M 


Uk 


4m  i — b  A tsmi  =  800  sec  =  13  min  20  sec, 

^  ’  Tj 

i=  1  1 


(10.1) 

(10.2) 

(10.3) 


respectively.  Thus  sine  dwell  testing  takes  significantly  longer  than  single  multisine  or  se¬ 
quential  multisine  testing.  However,  single  multisine  testing  cannot  be  used  because  of  the 
harmonic  distortion  that  results  from  system  nonlinearities. 

The  data  obtained  from  sine  dwell  excitation  is  used  to  calculate  the  total  harmonic 
distortion  of  DOT’s  dynamics  from  the  first  actuator  of  the  first  mirror  to  the  piston  po¬ 
sition  of  the  first  mirror.  The  total  harmonic  distortion  is  the  ratio  of  output  power  at  all 
super  harmonic  frequencies  to  the  output  power  at  the  base  frequency.  Figure  8  shows  that 
total  harmonic  distortion  of  DOT’s  dynamics  from  the  first  actuator  of  the  first  mirror  to  the 
piston  position  of  the  first  mirror.  Figure  8  demonstrates  that  there  is  significant  harmonic 
distortion  up  to  10  Hz  and  lower  levels  of  distortion  at  higher  frequencies. 

Figure  9  and  10  show  the  magnitude  and  the  phase  of  the  frequency  response  data 
obtained  from  the  single  multisine  and  sequential  multisine  excitation  signals.  The  two 
frequency  response  differ  most  significantly  between  0.5  and  30  Hz.  This  difference  can  be 
attributed  to  the  harmonic  distortion  that  is  present  in  the  single  multisine  data  but  is  not 
present  in  the  sequential  multisine  data. 
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Figure  8:  Total  harmonic  distortion  from  the  control  on  the  first  primary  mirror  segment  to 
the  piston  response  of  the  first  primary  mirror  segment. 


Figure  9:  The  magnitude  of  the  frequency  response  data  obtain  using  the  sequential  multisine 
excitation  (solid)  and  the  single  multisine  excitation  (dashed). 
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Figure  10:  The  phase  of  the  frequency  response  data  obtain  using  the  sequential  multisine 
excitation  (solid)  and  the  single  multisine  excitation  (dashed). 

11  Conclusions 

Many  excitation  signals  used  to  collect  frequency  response  data  for  linear  system  identifi¬ 
cation  are  poorly  suited  to  systems  that  exhibit  nonlinear  responses.  Specifically,  broadband 
excitation  signals  can  create  undesired  nonlinear  distortions  in  the  frequency  response  data 
that  are  generally  indistinguishable  from  the  linear  frequency  response.  Sine  dwell  excitation 
signals  allow  the  base  tone  response  to  be  separated  from  higher  order  harmonics.  However, 
sine  dwell  testing  can  be  very  time  consuming.  In  this  paper,  we  presented  a  novel  sequence 
of  multisine  excitation  signals  that  are  more  time  efficient  than  sine  dwell  excitation  while 
eliminating  the  effects  of  nonlinearity  induced  harmonics  from  the  frequency  response  data. 
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